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FOREWORD 


The  present  study  is  part  of  a  program  of  "Theoretical  Research  on 
Combustion  and  Wave  Propagation  in  Metal  Dust-Oxidant  Mixture^  being  con 
ducted  by  the  Division  of  Engineering,  Brown  University,  under  United 
States  Army  Ordnance  Contract  No.  DA-19-020-0RD-4761.  The  work  was  admin 
istered  by  the  Pyrotechnics  Laboratory  of  the  Picattiny  Arsenal,  Dover, 


N.  J. 


SUMMARY 


The  thermodynamic  properties  of  a  loosely-packed  granular  aggregate 
are  calculated.  An  application  of  these  formulas  to  the  study  of  a  shock 
wave  in  a  granular  mixture  is  indicated. 


THEEMODYKAMICS  OF  A  LOOSELY-PACKED  GRANULAR  AGGREGATE 


1.  Introduction 

We  shall  be  concerned  with  a  loosely-packed  granular  aggregate,  the 
space  between  the  granules  being  filled  with  a  perfect  gas*  We  assume  that 
the  diameter  of  the  Individual  granules  is  much  smaller  than  the  smallest 
macroscopic  scale  of  interest  and  that  in  every  macroscopic  volume  there 
are  a  large  number  of  granules.  Under  such  conditions  the  granular  aggre¬ 
gate  may  be  considered  as  a  continuum  from  the  macroscopic  viewpoint. 

The  mechanical  properties  of  such  a  continuum jof  course,  depends  on 
the  microscopic  state  of  the  granular  aggregate.  When  the  average  distance 
between  the  neighboring  granules  is  large  compared  to  the  average  diameter 
of  the  granules,  the  medium  behaves  (from  the  macroscopic  viewpoint)  as  a 
fluid.  It  is  in  fact  a  "dusty  gas".  The  thermodynamics  behavior  of  a 
dusty  gas  has  been  discussed  in  ref,  1, 

Let  us  compress  the  medium  (for  example,  by  a  piston)  until  the 
granules  begin  to  pack  together.  Further  compression  will  be  resisted  not 
only  by  the  gas  pressure  but  also  by  the  granules.  Since  the  elastic 
moduli  of  the  granules  are  much  higher  than  the  elastic  modulus  of  a  gas, 
we  may  consider  the  granxlLes  as  being  rigid.  Part  of  the  energy  expended 
in  compressing  the  gas  is  stored  up  in  the  gas,  the  remaining  part  is  ex¬ 
pended  against  the  friction  between  the  granules  and  cannot  be  recovered. 

If  the  piston  is  now  retracted  by  a  very  slight  amount,  the  granules  will 
no  longer  participate  in  exerting  any  force  on  the  piston.  The  force  on 
the  piston  required  to  maintain  this  small  change  drops  sharply  to  a:  value 
equal  to  that  contributed  by  the  gas  pressure.  If  we  now  recompress  the 
medium,  the  granules  will  again  participate  in  resisting  the  piston  and 
the  force  required  for  this  compression  rises  sharply.  Further  compression 
will  ultimately  cause  the  granules  to  lock  into  one  another,  the  medium 
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then  behaves  as  a  porous  rigid  solid.  The  thermodynamics  of  a  porous  rigid 
solid  may  he  deduced  from  that  of  a  porous  elastic  solid  which  has  been 
treated  in  ref.  2. 

In  this  report,  we  shall  examine  the  thermodynamics  of  a  loosely- 
packed  granular  aggregate  of  rigid  granules  in  the  intermediary  range  be¬ 
tween  a  "dusty  gas"  and  a  "porous  rigid  solid",  a  range  in  which  the 
mediiim  exhibits  what  is  usually  called  "plastic  deformation". 


2.  Entropy  Production  in  a  Deformation 

We  have  seen  that  when  our  medium  is  compressed,  part  of  the  work  done 
on  the  medium  cannot  be  recovered.  Of  course,  this  is  because  work  must 
be  expended  in  overcoming  intermole cular  friction  and  is  transformed  into 
heat,  resulting  in  a  production  of  entropy  inside  the  medium.  In  the 
present  case  the  entropy  produced  can  be  easily  calculated. 

We  begin  by  considering  a  cylinder  containing  a  unit  mass  of  the 
gross  medium.  At  one  end  of  the  cylinder  is  a  movable  piston.  Let  us  de¬ 
note  the  area  of  this  piston  by  Ap  >  the  loading  on  it  by  (T'  (in  force 
per  unit  area)  and  the  height  of  the  piston  'bY  ,  Part  of  the  loading 
is  carried  by  the  granules j  we  shall  denote  this  part  by  O' '  ,  The  re¬ 
maining  part,  denoted  by  0“"  ,  will  be  resisted  by  the  pressure  in  the 
gas .  Thus 


CT"  ■=  CT' '  4  (j- " 


(2.1) 


This  may  also  be  illustrated  in  a  stress-strain  diagram.  (See  Fig.  2). 

Let  0^  and  be  the  values  of  O'  and  70  when  the  medium  Just  ceases 

to  become  a  "dusty  gas"  at  the  temperature  ,  We  introduce  a  mealtil^  of 
compressive  strain  defined  by 


€  * 


In  figure  2,  we  plot  vs 

represents  the  loading  curve. 


(2.2) 

L®  :  ;iin  an  Isothermal  compression.  ABC 
BDA  represents  a  typical  unloading  curve. 
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The  ordinate  FB  represents  the  total  load  O'  ;  the  ordinate  FD  represents 
the  part  cr'"  resisted  hy  the  gas  while  DB  represents  the  part  O' '  . 

The  work  done  in  compressing  the  gas  is  or  -  O' cl  V  where  clV  is 

the  increase  in  volume  of  the  cylinder.  Thus,  a  positive  work  is  expended 
when  the  volume  decreases.  The  work  stored  in  the  gas  is  -  and 

that  transformed  into  heat  is  —  CJ^*<JiV  .  Hence  the  entropy  produced  in  a 
compression  per  unit  mass  of  the  gross  material  is 


G-'dV 

T 


(2.3) 


where  dV  is  clearly  the  specific  volume  of  the  material.  Note  that  this 
formula  applies  only  if  dV  <  O  ,  since  the  unloading  process  is  reversible. 
An  estimate  of  the  entropy  produced  in  a  deformation  can  be  easily 
made  if  the  work  done  in  compressing  the  material  is  given  and  the  compo¬ 
sition,  density  and  temperature  of  the  mixture  are  known.  The  composition 
of  the  material  may  he  conveniently  described  by  the  mass  fraction  Y'  of 
the  granules  in  the  gross  material.  Or  it  may  be  described  in  terms  of 
the  mass  fraction  Y''  of  gas  in  the  gross  material.  Of  course 

Y'  Y"  =  I  (2.4) 


Furthermore,  if  the  density  of  the  granules  is  #  which  is  assumed  known, 

the  total  volume  occupied  by  the  granules  in  the  mixture  is  .  The 

total  volume  occupied  by  gas  in  the  mixture  is  then  V  -  Y/f '  conse¬ 
quently,  by  the  gas  law. 


O-'*  -  'RT /CV-  Y/'^') 


(2.5) 


where  R  is  the  gas  constant  and  T  is  the  temperature  of  the  medium.  The 
entropy  produced  in  a  deformation  is: 


cls,  = 


ciw 


T  V_  Y’/5>' 

where  dw  is  the  work  done  in  compressing  the  material. 


(2.6) 
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3.  A  Fundamental  Relation 


The  first  law  of  thermodynamics  states  that 

dV  =  du  (3.1) 

where  ci(^  is  the  heat  added  to  the  material  and  du  is  the  internal 
energy  per  unit  mass.  If  denotes  the  entropy  added  to  the  material 

as  a  result  of  heat  transfer  to  the  medium. 


clq^  =  t4S^  (3.2) 

But 

dw  -  -T-civ  -  -  C-"<]v  + 

(3.3) 


therefore: 


=  ciO  +  ^"4V 

Now  dS  =  ,  therefore 

TdS>  =  dU  t 


(3A) 


(3.5) 


Introducing  Helmholtz  force  energy  per  unit  mass  of  the  gross  material, 


A  =  U  -  TS 

we  have 

c^A  =  -  S'^T  - 


Consequently, 


G-*''  = 


QV 


s 

u 


A  - 


T 


bT 


4,  Internal  Energy  and  Entropy  of  The  Material 
From  (2,5)  and  (3*8)  we  have 

V  -  r'/s' 


(3.6) 

(3.7) 

(3.8) 

(3.9) 

(3.10) 


(^.1) 
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since  YV5>'  is  the  volume  occupied  by  the  granules  in  a  unit  mass  of  the 
gross  material,  it  is  convenient  to  introduce  the  notation 


=  YVf'  (“t-s) 

Integrating  (4.1),  we  obtain 

^  ("w 

V.- 

where  ^  ^  and  Vo  are  the  values  of  A  and  V  when  T-T,  and  X  ■=■  , 

(cf.  eq..  2.2).  cI^Lt'^  in  (4.3)  is  an  arbitrary  function  satisfying  the  con¬ 
dition  =  o  .  Substituting  (4.3)  into  (3*9)  and  (3.10)  we  find: 


and 

(^.5) 

Thus,  the  internal  energy  of  our  material  is  a  function  of  T  only.  This 
is  as  it  should  be,  because 


u  =  yV  +  Y"^" 


(4.6) 


'r>a.i 

/  O  .!( 

and  the  Internal  energy  IJ  of  the  granules  and  U  of  the  gas  are  both 

A 

function  of  Y  ,  olf-T)  can  now  be  calculated  from  (4.5)  and  (4.6). 
Thus, 


=  A. 


(‘*•7) 


5.  Shock  Wave  in  A  Loosely-Packed  Granular  Medium 

As  a  simple  application,  the  forgoing  thermodynamic  formulas  may  be 
applied  to  calculate  the  discontinuous  change  of  state  across  a  shock  wave 
in  a  loosely-packed  granular  medium.  Of  course,  we  shall  also  need  an 
expression  for  the  variation  of  the  stress  component  0“'  with  strain  6  • 

A  relation  between  0-/  and  €  can,  in  principle,  be  determined 
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experimentally  ly  carrying  out  a  compression  test  of  the  material  in  an 
evacuated  container.  In  general,  a  relation  of  the  type 

O-'  ^  FCc->  (5^) 

can  be  expected  where  FCfe)  is  independent  of  the  temperature  ,  provided 
that  T  is  much  smaller  than  the  melting  point  of  the  granules.  Further- 
more,  we  assume  that  T «  o  and  fffe)  -^oo  as  ^ &  ,  at  which  point 

the  granules  are  so  tightly-packed  that  no  further  deformation  is  possible 
without  cracking  them.  A$  a  rough  approximate  a  O''  \js  &  curve  of  the 
form  shown  in  Fig.  3  niay  be  used. 

The  basic  equations  governing  the  jump  across  a  shock  wave  are  then: 


=  2-.  +  (5.3) 

^  r*- 

U  ^  UtT)  (5*5) 

Ch-  =  0^'+  0^"  (5.6) 

=  F  C  I  -  fc/f)  (5*^) 

=  'RT/(V-V^)  (5.8) 


where  subscripts  "1"  and  "2"  refer  respectively  to  the  conditions  ahead  of 
and  behind  the  shock,  ^  denotes  the  density  of  the  gross  material  and 
U.  denotes  the  velocity.  Eqs,  (5.2),  (5*3)  and  (5.^)  are  respectively  the 
continuity,  momentum  and  energy  equations;  eqs.  (5.5)  and  (5.8)  are  the 
equations  of  state,  while  (5.7)  describes  the  irreversible  part  of  the 
stress -strain  relation. 

If  we  substitute  eqs.  (5.5)>  (5.6),  (5.7)»  (5.8)  into  (5*2),  (5.3)> 
(5.^)>  we  obtain  a  system  of  three  equations  involving  and 
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Solving  these  in  terms  of  j  'V'l  >  the  dynamic  state  of  the 

material  behind  the  shock  is  completely  determined.  In  particular,  CT^ 
can  be  computed  from  (5.6). 
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